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In recent years, reallistic unquenched QCD simulations have been carried out with various lattice actions. In this 
report, I explain the progress in theory and algorithms and some of the physics results. 



1. INTRODUCTION 



The hadronic states in QCD have various energy scales: the smallest scales are the light pseudoscalar meson masses 
and uik, and the medium energy scales other light hadronic masses of order Aqcdj and the larger energy scales 
are heavy hadron masses mu and ms- These mass scales ranges from 100 MeV to 5 GeV, giving a hierarchy of order 
50. If one wants to fully describe such a system with completely controlled systematic errors, one should ideally carry 
out lattice simulations 1) with sufficiently large volume V = L'^ and small pseudoscalar meson masses mps to cover 
the chiral dynamics of the tt meson and 2) with sufficiently large heavy quark mass mg and large cutoff scale in 
order to cover the heavy meson mass scale. However, in reality, due to the limitation of the present computer power 
our lattice simulation can marginally cover the energy scales as 



Even if one cannot totally cover the energy scale, one can extrapolate the simulation result to the physical result 
with the help of the chiral perturbation theory and the heavy meson effective theory. Therefore, the real practical 
problem is to achieve lattice computations with larger volume, smaller light quark mass, and finer lattice as much as 
possible so that the kinematic parameters lie with the range where the effective theory descriptions are valid. 

In addition, in order to obtain the hadronic matrix element of the operator O with high precision, we need to 
determine the renormalization factor Zo very accurately since the matrix element in the continuum theory can be 
related to that on the lattice as 



In recent years, there have been several new developments both in theoretical formulations and algorithms. With 
these developments, recent lattice simulations can cover much wider energy scale so that the extrapolation of the 
quark masses to their physical values can be carried out reliably with the help of the low energy effective theories. 

1.1. Lattice fermion actions 

Let us now comment on various lattice fermion actions. The staggered fermion or improved staggered fermions 
(AsqTad, HISQ) have been extensively used for large scale unquenched simulations. The biggest advantages are the 
small simulation cost due to the small degrees of freedom and the exact partial chiral symmetry at finite lattice 
spacing. Since it has extra 'taste' degrees of freedom, square or quartic rooted trick are needed. The Wilson fermion 
or 0(a)-improved Wilson fermion (=Clover fermion) are now widely used in recent unquenched simulations. The 
biggest advantage is that this action is theoretically most simple and the computational cost is reasonable. Although 
the chiral symmetry is violated at finite lattice spacing, it can be recovered in the continuum limit. The Domain- Wall 
fermion [l| is a Wilson fermion in 5 dimensions with a mass term which has kinks in the 5-th dimensions. The chiral 
symmetry violations are exponentially suppressed as a function of the lattice extent in the 5-th dimension N^. It is 
-/V5 times more costly than the conventional Wilson fermion. The overlap 0] fermion has an exact chiral symmetry 
Q on the lattice. This action is very costly but still feasible by the new method to fix topology during the hybrid 



L ^ < rriT^ < mps < nix < Aqcd < rric < niQ < a ^ < nib. 



(1) 
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Table I: Large scale simulation projects 



Group 


Fermion Action 


Uf 


a(fm) 


L(fm) 


(MeV) 


MILC [10] 


Asqtad 


2+1 


0.06, 0.09, 0.12,0.15 


2.4-3.4 


> 240 


PACS-CS [11] 


Clover 


2+1 


0.07, 0.09, 0.10 


3 


> 150 


BMW [13] 


Clover 


2+1 


0.065,0.085,0.125 


2-4 


> 300 


CERN [14. 151 


Wilson 


2 


0.052,0.072,0.078 


2.5-3.5 


> 300 


CLS/ Alpha [16] 


Clover 


2 


0.04, 0.06, 0.08, 


1.9-2.4 


> 300 


ETMC [12] 


tmWilson 


2 


0.07, 0.09, 0.10 


3 


> 300 


RBC/UKQCD [IT] 


Domain-wall 


2+1 


0.086,0.12 


3 


> 330 


JLQCD [18J 


Overlap 


2 


0.12 


2 


> 300 






2+1 


0.11 


2 


> 300 



Monte Carlo simulation. The twisted mass Wilson (tmWilson) is the Wilson fermion with chirally rotated mass. It 
has a partial exact chiral symmetry but partial vector symmetry is violated, instead. At maximal twist, this action 
is automatically 0(a)-improved. 



1.2. New algorithms in unquenched simulations 



It is known that the simulation cost grows towards smaller quark mass, larger volume, and finer lattice spacing. In 
addition to the obvious effect from simply having more lattice points, there are effects from more iterations in Dirac 
matrix inversion, smaller step size in keeping the acceptance rate constant, and larger Monte Carlo trajectories to 
compensate longer autocorrelations. For example, the empirical formula [J] for Nj = 2 QCD with 0(a)-improved 
Wilson fermion based on the studies by CP-PACS and JLQCD collaborations with conventional algorithm is 



cost [Tflops • years] — C 



T^conf 




0.6 


6 


L 


5 


"O.lfm" 


1000 








3 fm 




a 



(3) 



with C ~ 2.8 . From this formula, one can estimate that a typical simulation of 1000 configurations with a = O.lfm, 
L = 3 fm would cost 25 Tflops years for = 300 MeV and more than 600 Tflops years for physical pion mass. 
For this reason, it has been thought that the staggered fermion is the only feasible choice for unquenched QCD 
simulations for ~ 300 MeV with Tflops machine. 

One of the most striking developments in lattice QCD is the proposal of new algorithms to reduce the simulation 
cost of unquenched QCD simulation. The new algorithms are the Hybrid Monte Carlo based on the combination 
of the 'preconditionning' (mass preconditionning p, Q or domain decomposition 0, S]) to reduce the cost of the 
Dirac operator inversion and 'multi-time scale' in Molecular dynamics to reduce the frequency of the Dirac operator 
inversion in time-step. It was found that these new algorithms can give rise to significant cost reductions of factor 
20-30. Therefore, the study of 0(a)-improved Wilson fermions has also become quite feasible with 0(1-10) Tfiops 
machines. The same algorithm also enables large scale simulations of Domain-wall or overlap fermions, although one 
must restrict oneself to work on a single relatively coarse lattice spacing with the present computer. 



1.3. Nonperturbative renormalization 

Another important theoretical methods is the nonperturbative renormalization. There are several different methods. 
The first method is a simple method using the PCAC relation A^A^ = 2mqP. If the axial current is exactly 
conserved on the lattice (ex. staggered, tmWilson, Overlap), one can use PCAC relation to obtain the decay 

constant of the pseudoscalar meson with nonperturbative accuracy as /ps = lim — 1 — ^— . The second method is 

a^O "Tips 

the nonperturbative renormalization by the Regularization Independent Momentum (RI-MOM) scheme proposed by 
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Figure 1: Low energy constants from NLO SU(2) ChPT fits of the lattice data from various groups. For comparison, 
phenomenological estimates by Gasser and Leutwyler are also plotted. 

Table II: SU(3) breaking in the pseudoscalar decay constant. 



Group 



Action 



input fR-ffT, 



MILC 2+1 Asqtad i^r,7v,K 

HPQCD/MILC 2+1 Asqtad+HISQ mr,^,K 

PACS-CS 2+1 Clover 

BMW 2+1 Clover 

RBC/UKQCD 2+1 Domain- waU 

JLQCD 2+1 Overlap 

ETMC 2 tmWilson 



1.202(3)(: 
1.189(7) 
1.19(2) 
1.19(1)(1) 
1.22(2)(6) 
U,m^,K 1.20(3) 
U,m^,K 1.196(13)(7)(8) 



ma.Tv,K 



Martinelli at al. [l^. This is a scheme defined by the quark-gluon amplitudes in Landau gauge with large off-shell 
momentum in Euclidean region. Schrodinger functional (SF) scheme proposed by the Alpha collaboration 
This scheme is defined the amplitudes in a box with boundaries with physical box size under the external gauge 
potential introduced by the boundary values of the gauge field. 



2. LIGHT FLAVOR PHYSICS 



2.1. Quark mass dependence of the pseudoscalar masses and decay constants 

Next-to-Leading Order (NLO) SU(2) Chiral Perturbation Theory (ChPT) predicts the quark mass dependences 
of the pseudoscalar masses and decay constants as 



'PS 



fps — f 



1 



2m,gB 



u 



(4) 



where raT^= 139.6 MeV. The ratio of the decay constant is also obtained by various groups. Table |TT] lists the 
ratio Jk/ f-K in 2 and 2+1 flavor QCD from various groups. A remarkable observation was made by the 2+1 flavor 
simulation PACS-CS collaboration llj, where the lightest quark mass is as small as the physical up and down quark 
masses. They find that the non-analytic behavior as functions of mud arising from the chiral log can be clearly 
reproduced for m^/mud and fx / fix- RBC/UKQCD collaboration found that the NLO partially quenched chiral 
perturbation theory (PQChPT) with the expansion parameter x — m^/(47r/)^ (a:-expansion)fits the data only for 
< 450 MeV. While JLQCD collaboration finds similar results for ChPT, they also find that the NLO ChPT 
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Figure 2: ChPT fit with x-expansion and ^-expansion, x-expansion makes the convergence of the ChPT fits better. 

with expansion parameter ^ = / (AirfT^)'^ (^-expansion) makes the convergence better for heavier quark masses and 
NNLO ChPT with ^-expansion can nicely describe the lattice data up to the strange quark mass regime as shown in 
Fig. [21 [18]. 



2.2. B 



K 



Indirect CP violation in the K meson system eK is one of the most crucial quantities to test the standard model 
and the physics beyond. The experimental value is determined with high accuracy as 



lexl = (2.233 ±0.015) x 10 



-3 



(5) 



Theoretically this quantity is described as |e/^| = f[p,'q) x C{^i) x Bk{^))- Here, f[p,'q) is a factor which depends 
on the CKM matrix elements, C(/i) is the Wilson coefficient from short-distance QCD corrections and Bk{ij.) is the 
bag parameter defined as 



{K°\ [drji - i5)sdj,{i - 75)^] (/i)i^°) 

8 f2 2 



(6) 



The main problem in unquenched lattice calculations is the mixing of operators with wrong chiralities or tastes. For 
Wilson-type fermions there exists mixing with operators of the same dimension having wrong chirality. Since they 
cannot be extrapolated away in the continuum limit, the only way to remove the contamination is to determine the 
counter term nonperturbatively. However, due to the chiral enhancement of the operators with wrong chirality, it is 
quite difficult to control the systematic error. The operator mixing for staggered fermion is also quite complicated 
and only perturbative renormalization exists. The overlap fermion is free from operator mixing owing to the exact 
chiral symmetry, while the operator mixing for the domain-wall fermion is exponentially suppressed and practically 
under control. In fact, for both cases studies with nonperturbative renormalization in RI-MOM scheme, one cannot 
see visible effect of the the operator mixing. For tmWilson, due to the parity odd nature and the cs symmetry 
the operator mixing can be prohibited for 2-flavor QCD. Therefore, tmWilson is another promising approach to the 
precise determination of Bk- However, for 2-1-1 flavors it is difficult to full realize the 0(a)-improvent, no mixing 
and non-degenerate (771^,771^) quark simultaneously. 

HPQCD collaboration studied Bk on 2+1 flavor MILC configuration with Aqstad (staggered ) sea quark with 
lattice spacing a = 0.125 fm 23[. They employ two different improved staggered actions (Asqtad and HYP) for the 
valence quark. The calculation was carried out by degenerate valence quarks for the K meson {mi, TO2) = rasjl while 
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the sea quark is chosen to be 0.2msea, 0.4ms. Chiral extrapolation is performed by linear fit. The renormalization 
of AS = 2 operator is obtained by perturbation at 1-loop. They find B = 0.83(18), where the dominant error comes 
from the unknown higher order perturbative correction of the renormalization factor of 0(0^) ~ 20%. RBC/UKQCD 



collaborations computed Bk with domain- wall fermion in 2+1 flavor QCD at lattice spacing a — 0.11 fm [17[. They 
have 3 points for the sea quark and 7 combinations of valence quark masses (mi, 777.2). They fit the data with 
next-to-leading order partially quenched chiral perturbation theory (NLO PQChPT). The renormalization factor 
is determined nonperturbatively by RI-MOM scheme. They obtain B = 0.720(13)stat.(17)sys., where the dominant 
error comes from the discretization error of O(a^) — 4%. JLQCD collaborations study Bk with overlap fermion 



in 2 flavor QCD at lattice spacing a = 0.12 fm on physical volume with L — 2/m [2J]. They have 4 points for 
the sea quark and 10 combinations of valence quark masses (mi,m2). They fit the data with NLO PQChPT. The 
renormalization factor is determined nonperturbatively by RI-MOM scheme. They obtain B = 0.734(5)stat.(50)sys., 
where the dominant error comes from the finite size effect of order 5%. 

It should be noted that the long standing operator mixing problem is solved with the advent of Ginsparg- Wilson 
fermions (domain-wall, overlap) with (almost) exact chiral symmetry. Thus the above studies are the real start of 
the precision study of Bk for which significant progress will be expected near future. RBC/UKQCD collaboration 
are planning to work at finer lattice spacing a = 0.09 fm to reduce the discretization error. JLQCD collaborations 
are studying 2+1 fiavor QCD for physical volumes of 2 fm and 3 fm. ETMC and other groups are studying on mixed 
action approach where they employ the Ginsparg- Wilson valence quarks on gauge configuration with tni Wilson or 
staggered dynamical quarks. This can also be an economical and promising approach, with the help of mixed action 
PQChPT so that one can even take the continuum limit with the existing 2+1 fiavor gauge configurations. 

More precise determination of Bk are important in view of the recent tension in sin(20i) versus the determination 
of unitarity triangle by ex and AuibJAb^ as pointed out by Lunghi and Soni ^] or Buras and Guadagnoli [2^. 



2.3. Nucleon sigma term 

Nucleon sigma term is the finite quark mass effect of the nucleon mass and is defined by the following matrix 
element 

cr^^ = ■mud{N\uu + dd\N) , (7) 

which contains both the valence and the sea quark contributions. The sea quark contribution also arise from the 
strange quark and it is often parameterized by the ratio y defined as 

^ {N\uu + dd\N) ^ ' 

The sigma term and the strange quark content can be determined using the quark mass dependence of the nucleon 
mass. The Feynman-Hellman theorem tells that the quark mass dependence can be related to the sigma term as 

= (7V|7277 + Jd|iV)eo„n, ^ = (iV| 7/77 + Jd| TV) disc , (9) 

where "conn" and "disc" denotes contributions to the sigma terms. The disconnected contribution is nothing but 
2{N\ss\N) when the sea quark mass is equal to strange quark mass. 

JLQCD collaboration studied the quark mass dependence of the nucleon mass in 2-flavor QCD simulation with 
dynamical overlap fermion [27]. Fitting the nucleon mass for the combination of 6 sea quark masses and 9 valence 
quark masses by partially quenched baryon chiral perturbation theory. They found that the by fitting the data with 
BChPT in finite volume finite size effect for the baryon mass in (2/m)'* box is significant. After correcting the finite 
size effect, they obtain 

a^N = h2{2){tf){tl) MeV. (10) 
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Figure 3: Left panel shows the quark mass dependence of the nucleon mass in 2-flavor QCD. Right panel shows the connected 
and disconnected contributions to the sigma term. 



Fitting the nucleon mass data with PQChPT and differentiating the nucleon mass with valence and sea quark 
masses. They found that the disconnected contribution is much smaller than connected one. As a result the semi- 
quenched estimate of the y-parameter is 

y ^ QmG{l&)tl)t\). (11) 

Previous lattice results y-parameters using Wilson-type fermions are 0.66(15), 0.36(3) for quenched QCD and 
0.59(13) for 2-flavor QCD, which are significantly larger than the above result. It was pointed out that the the additive 
quark mass shift with Wilson-type fermion can give a significant contamination to the disconnected contribution. 
After subtracting this artifact, UKQCD obtained y = —0.28(33) in 2-flavor QCD. Since they subtract the power- 
divergent contribution numerically, the result suffer from large statistical error. In contrast, the study with overlap 
fermion by JLQCD do not suffer such artifact owing to the exact chiral symmetry. 

The exact chiral symmetry on the lattice enables novel studies of chiral dynamics. It can also open up new directions 
which has not been possible in other approaches. Examples of such studies are given by JLQCD collaboration [isj. 



29, 30, 31, 3 



not been possi 



3. CHARM QUARK PHYSICS 
3.1. D meson decay constants 

Since the charm quark mass is of order 1 GeV, the discretization errors are larger than those for up, down, and 
strange quark if one take the same approach used for the light quark. In order to make a reliable extrapolation 
should make the lattice one has to employ either finer lattices or improved action. 

The first choice is made by the Alpha collaboration. They carry out a quenched study the Ds meson decay 
constant on lattices with lattice spacings a — 0.04,0.06,0.08,0.1 fm, which correspond to the discretization errors of 
[anicY — 0.04,0.09,0.16,0.25 36]. Using the 0(a)-improved Wilson fermion and nonperturbative rcnormalization 
in SF-scheme, they obtain 

flf=^ = 252(2) MeV. (12) 

They same strategy is in progress for = 2 QCD by the Coordinated Lattice Simulation(CLS)/Alpha collaboration. 
This is a quite promising approach but requires dedicated effort. 

Improved fermion actions are also effective in reducing the discretization error. One of the popular approaches is 
the fermilab formalism and other related formalism in which one make a mass-dependent re-interpretation of 0{a)- 
improved Wilson fermion. While it is possible to achieve non-perturbative improvement, mass-dependent errors of 
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the leading operators are removed at tree and 1-loop in practice. The dominant short distance error A for the charm 
quark system is 

A = Coal + cias{aAQci))iamc) + C2(aAQCD)^ H (13) 



Highly Improved Staggered Quark (HISQ) is another improved action [37|. Since it is based on the standard 
Symanzik improvement which treats the discretization errors in powers of (ap)^ and (awg)^, standard nonpertur- 
bative renormalization technique also applies. The advantage of this approach is that the discretization errors are 
restricted due to the exact chiral symmetry of the staggered fermion. In general, after carrying out nonperturbative 
renormalization the dominant short distance error A for the charm quark system can be expanded in powers of 
(oTOc)^. In HISQ action, they remove the leading 0{a^) error and the taste changing interaction at tree level so that 
the dominant short distance then becomes 

A = ci(amc)'* + C2as(amc)^ H . (14) 

In principle, HISQ has nonperturbative accuracy in the continuum limit. Whether typical lattice spacings a ~ 0.10 
fm are small enough to have a good scaling for charm quark system as expected from naive order counting should 
be explicitly examined. 

HPQCD collaboration studied charm quark system on MILC lattices with 2+1 dynamical Asqtad quark using 
HISQ action for the valence quark with the lattice spacings a = 0.09,0.12,0.15 fm [38|. The chiral and continuum 
extrapolation for the pseudoscalar and masses and decay constants is carried out by the global fit of the data using the 
ChPT and heavy meson effective theory formula including ©(a^a^) and 0{a'^) errors. The scale is set by the upsilon 
spectrum and m„d, f^s and rric are determined from m^, mx and m^^. Since there are no more free parameter, 
while fjj^ and fn^ are pure predictions of theory, it should be noted that f^, fx, tid^, m^i^ and ttt-d^ can also 
give non-trivial consistency check. It turns out mixed action approach with HISQ valence and Asqstad sea quark 
can control the discretization error in the charm quark system nicely. They showed that the D and Dg meson 
masses are reproduced in the continuum limit. The pseudoscalar decay constant are computed with nonperturbative 
renormalization by AWTI, the decay constants /^r and fx agree with experiment very well. Using the same technique 
they obtain fu^ and fo, with 1-2% error which is dominated by the continuum and chiral extrapolation. 

fo, = 207(4) MeV, fo^ - 241(3) MeV, /dJId, = 1.164(11). (15) 

Fermilab/MILC collaboration also studied the D meson decay constants on the same 2+1 flavor MILC configura- 
tions 3^ 49 1 . They used AsqTad action for the light valence quark and fermilab action for the charm quark. The 
data are fitted by the PQstaggered ChPT which include full NLO and analytic part of the NNLO contributions and 
O(a^) errors. The heavy- light axial vector current is renormalized by partially nonperturbative renormalization in 
which they use the renormalization factors of the heavy- heavy and light-light vector currents Zy^ , Zyj^ and evaluate 
the remaining perturbative correction p^"^ 



They expect that vector current renormalization factors capture the dominant contribution from the wavefunction 
renormalization. Their preliminary result is 

fn, = 207(11) MeV, fo^ = 249(11) MeV, IdJId, - 1-20(3). (17) 

ETMC studied the D meson decay constants for Uf ~ 2 QCD using dynamical tmWilson with lattice spacing 
a = 0.67, 0.86, 0.10 fm 40]. Since they work at maximal twist which guarantees automatic 0(a)— improvement, the 
leading discretization error is O(a^). They take degenerate ud quark masses for the sea and valence quark in the range 
of ruud = (0.2 — 0.4)ms. The data are fitted with SU(2) heavy meson ChPT with O(a^) term and nonperturbative 
renormalization by AWTI is used for the heavy- light axial vector current. They obtain the preliminary result as 

fD, = 197(16) MeV, fo^ = 244(12) MeV, fajfo, = 1.24(5). (18) 
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Fig. |4] shows the comparison of the D meson decay constants. It can be seen that f^a from HPQCD is in agreement 
error with experiment, while f^,^ has 3 a deviation from CLEO result. Two recent preliminary results from Fermi- 
lab/MILC and ETMC give consistent result. It was pointed out by Becirevic et al. (4l| that possible experiment 
error by misidentifying D decay could only give an effect smaller than a few percent. As was pointed out by 

the lattice talk by Fermilab group, the 3 a deviation is dictated by the error from CLEO even if we enlarge the error 
from HPQCD by factor 3, the deviation still 2.8 a. Therefore the fo puzzle still remains. Further theoretical and 
experimental studies are required. Combining HPQCD(n/ = 2 + 1), FNAL/MILC(n/ = 2 + 1), and ETMC(n/ = 2) 




O CLEG 

□ ALPHA 2003 

♦ ETMC 2008 

▲ HPQCD 2008 

■* FNAL/MILC 2008 

I 1 , , , U \ , , I I , , 

150 200 250 300 350 

f„.f„ (MeV) 

Figure 4: D meson decay constants 
results where I include additional 5% error for Uf = 2 results my world average is 

fo, = 206(4) MeV, /d, = 243(3) MeV, (19) 



N=0 
N=2 



N=2+l 



3.2. Charm quark mass 

Another important progress is the precise determination of the charm, quark mass. HPQCD and Chetyrkin, Kuhn, 
Steinhauser and Sturm 4^, |43| determined nic and as from the moments the pseudoscalar correlators 

^"^Ef-) ^W' G(i)^a6Es(amS)2(0|j5(x,t)j5(0,0)|0), js^^cTsV-- (20) 



t ^ 



They match the continuum limit of the lattice calculation of Gn with the continuum perturbation theory of Gn at 
4-loop. The lattice calculation was performed with HISQ valence charm based MILC lattice for ri/ = 2 + 1 QCD at 
lattice spacings a = 0.09, 0.12, 0.15 fm. 

In order to reduce the discretization error, modified ratios i?„ are defined 

i?4 = Gi/Gf\ Rn = 1^ (g„/g1°)) '^^"^'^ for n > 6, (21) 

where Gn'^ are the n-th moments on the lattice at the tree level. The lattice spacing dependence of the moments can 
be fitted nicely with quadratic functions in a^. Averaging for n — 6,8, 10 one can determine the charm quark mass 
and the strong coupling as 

mf^(3GeV") = 0.986(10), ^^(M^, ti/ = 5) = 0.1174(12). (22) 

The key ingredients are the O(a^) improvement and the exact chiral symmetry, which can control the discretization 
error and also allow nonperturbative renormalization by AWTI. 
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Table III: Quenched QCD results of /s and mj, 



Group 


Method 


/s, (MeV) 


mf ^(mt) (GeV) 


Alpha [44,45] 
Rome II [46, 47] 
Guazzini et al. [481 


Static + 1/mb 

finite volume step scaling 

combination 


193(7) 

192(6)(4) 

191(6) 


4.35(5) 

4.33(10) 

4.42(6) 



3.3. B PHYSICS 

Even with the finest lattice spacing for the present unquenched lattice calculation a = 0.04 fm, the bottom quark 
mass in lattice unit is in the range of am?, ~ I. Therefore, HISQ action completely looses control over the discretization 
error. Therefore, the only practically available methods to study B physics are either to use the effective action such 
as NRQCD or fermilab a action or to make a very precise computation in both the static limit and in the charm 
quark mass regime and make an interpolation. While the former approach are taken by HPQCD and FNAL/MILC 
collaborations and has already produced many results in the past few years, there has also been remarkable progress 
in the latter approach, where serious feasibility tests of the theoretical method have been carried out in quenched 
QCD. 

Alpha collaboration worked on the lattice HQET with nonperturbative accuracy 



44 



45l |. Their idea is to make a 



nonperturbative matching of the HQET and QCD at very fine lattices in small volume and then evolve the lattice 
HQET to coarse and larger lattice by step scaling. They could also work on l/rrih correction in HQET. 

Rome II group used step scaling using 0(a)-improved Wilson fermion with nonperturbative renormalization 4^ 47 1. 
They make nonperturbative calculations of B meson directly at bottom quark mass in a small volume L = 0.4 fm 
at very fine lattice spacing. Then evaluate the finite volume correction with coarser lattice by extrapolations from 
smaller quark mass regime. 

Guazzini, Sommer and Tantalo combined the above two methods and used the static results to evaluate the finite 
volume correction by interpolation 



The quenched QCD results are summarized in Table [1111 It is remarkable that one can control the systematic error 
from the bottom quark and achieve 2-3% accuracy in the matrix element, ny = 2 unquenched QCD studies are in 
progress by Alpha collaboration. 

FNAL/MILC collaboration recently update their = 2 + 1 calculation of fs^ and fs^ |49|- Fig. [5] shows the Bs 
meson decay constants in quenched QCD and = 2, 2 + 1 unquenched QCD simulations. Using the input such as 
TOp or To = 0.5 fm for rif =^ 0,2, one finds that /s^ increase from quenched QCD to ny = 2 + 1 in contrast to f^^ . The 
SU(3) breaking of the fs in 2+1 flavor QCD is obtained by HPQCD and FNAL/MILC collaboration also reported 
their result at Lattice 2008. Since FNAL/MILC results of fs are still preliminary, I just quote their numbers and 
do not take the world average. 



IbJIb, = 1.20(3) HPQCD, JeJIb, = 1.25(4) HPQCD. 



(23) 



3.4. B D*lv 

FNAL/MILC collaboration recently made an update of B — > D*lv form factor determination based on 2+1 MILC 
lattices at lattice spacings a = 0.09,0.12,0.15 fm [50]. They employ the Asqtad action for the light sea and valence 
quarks and the fermilab action for the heavy quark. They developed a new double ratio method to determined the 
B — > D*lv form factor as 

I ^ |2 _ {D*\nn^h\B){B\h,^,c\D*) 
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Figure 5: Bs meson decay constant. 



They carry out the chiral extrapolation and continuum extrapolation using the NLO ChPT and analytic NNLO 
staggered ChPT formula. They obtain 

J^"^=2+i(l) = 0.921(11)(19) (25) 

as a preliminary result. 

Rome II group studied B D*lv form factor in quenched QCD using finite size step scaling method and combined 
with the HQET calculation 51|, 53] ■ They used twisted boundary condition to get small continuous momentum 
recoil. As a result they could obtain the slope of the w dependence. They obtain 

jr«/=0(i) ^ 0.917(8)(5). (26) 

This proves that the finite size step scaling method is very promising also for the form factor calculation. Applications 
to unquenched QCD simulation are awaited. 



3.5. B nlu 

FNAL/MILC collaboration made a new analysis oi B —^ ttIv form factor using MILC configuration of 2+1 flavor 
QCD at a = 0.09,0.12,0.15 fm {Hsll . Parameterizing the the form factors as 

f^{E^)p^ = {7r\V,\B)/V2^, f\\{E^) = {n\Vo\B)/V2^ (27) 

and using NLO staggered ChPT 

fl± = Co [l + chiral logs + mi + C2(2m„ + to^) + cg-E^ + C4,EI + cga^] + 0{ml,El), (28) 

they carried out the continuum and chiral extrapolation. They then used the combined z fit of Babar and lattice 
data based on the dispersion relation formula and extracted \Vab\ as 

iKbl = 2.94(35) X 10-3. (29) 

3.6. New applications 

The nonperturbative renormalization developed for QCD simulation for precise test of the standard model is also 
useful for exploring the beyond standard model physics. Recently, the running coupling constant in Schrodinger 
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functional scheme is applied to QCD with many flavors. Appelquist et al. 5J| found an evidence of an infrared fixed 
point in strong coupling regime in QCD with twelve flavors. This may serve to construct an example for the walking 
technicolor. Application to more realistic models are now in progress. 



4. CONCLUSIONS 

In conclusion, the recent developments of unquenched lattice simulation allows us to study the realistic 2+1 flavor 
QCD. A lot of nonperturbative computation in light flavor physics is now well under control. Highly improved 
fermion actions are beginning to control the charm physics. The b quark still suffers from large systematic errors 
but several promising methods have been tested and giving precise determinations of weak matrix elements of the B 
mesons in quenched QCD. It was found that the exact chiral symmetry is quite useful in some cases. 

Coordinated work to combine advance techniques will lead to few percent accuracy within few years, including 
cross-checks with different lattice actions. 
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